Wigner formula of rotation matrices and quantum walks 
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Quantization of a random-walk model is performed by giving a qudit (a multi-component wave 
function) to a walker at site and by introducing a quantum coin, which is a matrix representation 
of a unitary transformation. In quantum walks, the qudit of walker is mixed according to the quan- 
tum coin at each time step, when the walker hops to other sites. As special cases of the quantum 
walks driven by high-dimensional quantum coins generally studied by Brun, Carteret, and Ambai- 
nis, we study the models obtained by choosing rotation as the unitary transformation, whose matrix 
representations determine quantum coins. We show that Wigner's {2j -\~ l)-dimensional unitary rep- 
resentations of rotations with half-integers j's are useful to analyze the probability laws of quantum 
walks. For any value of half-integer j, convergence of all moments of walker's pseudovelocity in the 
long-time limit is proved. It is generally shown for the present models that, if {2j -\- 1) is even, the 
probability measure of limit distribution is given by a superposition of {2j -\- l)/2 terms of scaled 
Konno's density functions, and if (2j -f 1) is odd, it is a superposition of j terms of scaled Konno's 
density functions and a Dirac's delta function at the origin. For the two-, three-, and four-component 
models, the probability densities of limit distributions are explicitly calculated and their dependence 
on the parameters of quantum coins and on the initial qudit of walker is completely determined. 
Comparison with computer simulation results is also shown. 

PACS numbers: 03.65.-w,05.40.-a,03.67.-a 
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I. INTRODUCTION 



No one doubts the importance of random-walk mod- 
els in physics, mathematics and computer sciences. In 
particular, when we explain basic concepts of statisti- 
cal physics [l[ , stochastic processes in physics and chem- 
istry and stochastic algorithms [3|, introduction of 
random-walk models is very useful and effective. It is 
interesting to see that systematic study of quantization 
of random walks is not old 

0, H, B Q- As expected, 
the study of quantum walks is fruitful and its results 
have been applied to solve the transport problems in 
solid-state physics of strongly correlated electron systems 
to derive non-Gaussian- type central limit theorems 
in probability theory [§, [l^, [Hi, and to invent new al- 
gorithms in the quantum information theory [l3| . 
The research field of quantum walks is growing widely 

[S E H [O, E [H 

and mathematical understanding 
of the new models is becoming deeper [lO, HH, H^, IH, 
in recent years. 



observation of physical quantities, 



when we calcu- 



It should be noted here that, from the view-point 
of standard quantum mechanics, "to quantize random 
walks" is a contradictory concept, since in quantum me- 
chanics, time-evolution of a state vector or a wave 
function ^'(x, f), is given by a deterministic unitary trans- 
formation associated with the Hamiltonian and probabil- 
ity concept appears in the theory only when we perform 
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late the probability density p{x, t) — \^{x, t)p at a given 
time. On the other hand, random walk is a typical ex- 
ample of stochastic processes, in which we toss a coin to 
select a walk at each time step. In an earlier paper ^2^, 
it was shown that the one-dimensional standard random 
walk can be realized by a random-turn model [26j . in 
which a coin is represented by a 2 x 2 stochastic matrix 
and that, if we replace the matrix by a 2 x 2 unitary ma- 
trix, a one-dimensional quantum- walk model is obtained. 
This argument is not only heuristic but also generic, since 
it implies that quantization of random-walk models can 
be done by introducing appropriate unitary matrices such 
that they play the roles of "quantum coins" . The ob- 
tained time-evolution of quantum walk is described by 
multi- component version of quantum-mechanical equa- 
tion of motion. For the standard quantum-walk model 
on one-dimensional lattice Z — {• • • , —2, —1, 0, 1, 2, • • •} 
with the nearest-neighbor hopping, the equation is iden- 
tified with the Weyl equation for two-component wave 
functions [1^. It should be remarked that such multi- 
component equations have been usually used in relativis- 
tic quantum mechanics [27| . 

It should be noted that in order to discuss the rela- 
tionship between classical and quantum walks Brun et 
al. showed a method to construct discrete quantum- 
walk models for multi-component wave functions driven 
by high-dimensional quantum coins that are greater than 
2x2 matrices There they considered a 2^- 

dimensional space with M — 2,3, ■ ■ ■, whose basis is given 
by tensor products of M binary-vectors. A quantum coin, 
which is first defined as an abstract unitary transforma- 
tion U, has a 2^ x 2^ matrix representation in the space 



[18|. Their method is very general and useful as well as 
the tensor-product method is in the group representation 
theory [ll]. See also Refs. [l3,[2l,[3^ for tensor-product 
models. 

In the present paper, we adopt the rotation operator R 
specified by the Euler angles a, f3, 7 (see Eq.Q in Sec. II) 
as the unitary transformation U and introduce a family 
of quantum- walk models on one-dimensional lattice Z. 
The two-dimensional representations of R can be identi- 
fied with 2x2 quantum-coin matrices used in the stan- 
dard quantum walks with qubits. In the wave-number 
space (fc-space) the shift matrix S{k) is generally given 
by a matrix representation of R, if we set the parameters 
a = -2k, (3 = 7 = (see Sec.III). Our choice U = R 
is thus very suitable to study one-dimensional quantum 
walk problems. 

Instead of using the tensor-product representations fol- 
lowing Brun et at, we will use in this paper the {2j + 1)- 
dimensional unitary representations (a, /3, 7) of the 
rotation operator R with half- integers j's as quantum 
coins, which are called the rotation matrices [3l|. The 
reason is the following. The 2*^-dimensional represen- 
tations with M > 2 constructed by tensor products are 
reducible; By changing basis through appropriate orthog- 
onal transformations, tensor-product matrices are block- 
diagonalized. Wigner's theory showed that each block 
with size 2j + 1 {j — 0, 1/2, 1, 3/2, • • •) is given by a rota- 
tion matrix R^^^^ {a, (3, 7) [32|- In other words, we will use 
the irreducible representations of R. In Sec. IV, reducibil- 
ity of the tensor-product models of quantum walks to the 
present models will be generally shown and demonstrated 
using examples. 

In our model with j and three parameters (a,/3,7), a 
quantum walker is assumed to be at the origin at time 
t = with a (2j + l)-component qudit 



( \ 



\ 1-0 I 



with 



3 

E 



(1) 



where G C (complex numbers), m — —j, ~j -I- 1, • • • , j. 
At each time step t = 1, 2, 3, • • •, the components of qudit 
are mixed according to a quantum coin R^^\a,l3,-i) and 
the walker hops to (2j -I- 1) sites on Z, as illustrated by 
Figfflfor j ^ 1/2,1,3/2 and 2 (see in Sec.III). 

When j = 1/2, R'^^/'^\a, can be identified with 
an element of SU(2) appropriately parameterized by the 
three variables (the Cayley-Klein parameters) and the 
model is reduced to the standard two-component model 
[2^ . It should be noted here that, when j is an integer 
(i.e. (2j -I- 1) is odd), the walker can stay at the same 
position in a step. 

Using the method of [ll|, we will prove that any mo- 
ments of pseudovelocity of walker, which is defined by 
Xt/t (the position at time i, Xt, divided by t), converges 
in the long-time limit t ^ 00, and show that the probabil- 
ity measure of limit distribution is generally described by 
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FIG. 1: (Color online) Elementary hopping of quantum 
walker in the models with (a) j = 1/2 (two-component 
model), (b) j = 1 (three-component model), (c) j = 3/2 (four- 
component model), and (d) j = 2 (five-component model). 
When (2j -I- 1) is odd, the walker can stay at the same posi- 
tion in a step, as shown by cyclic arrows at the origin in the 
cases (b) and (d). 



a superposition of appropriately scaled forms of a func- 
tion 



r(l — x'^)Va? 



-2 -fl^KI'^l}' 



(2) 



where i{ui} denotes the indicator function of a condition 



1 if w is satisfied and 1 



otherwise, 



and a is a real parameter. It is the density function 
first introduced by Konno to describe the limit distribu- 
tions of the standard two-component quantum walks in 
his weak limit-theorem [l,[l3]- (As shown in FiglUa) in 
Sec. VI, /i(x; a) is inversed bell-shaped on a finite support 
X e (—a, a) in big contrast with the Gaussian distribu- 
tion, which describes the diffusion scaling limit of the 
classical random walks.) More precisely speaking, when 
(2j -|- 1) is even, the probability density function of limit 
distribution consists of (2j -I- 1)/2 terms of Konno's den- 
sity functions and when (2j + 1) is odd, it does of j 
terms of Konno's and a point mass at the origin, which 
corresponds to the positive probability to retain the po- 
sition of walker in a step (see Eg. ((72)) ). 

The weight functions Al^-''™) of these Konno's density 
functions (and the weight A^-') of Dirac's delta function at 
the origin, when (2j -I- 1) is odd) in the superposition de- 
pend on the parameters of quantum coin and the (2j -I- 1) 
components of initial qudit ([T]) of quantum walker. We 
have found that the representation theory of groups [2^ 
is useful to calculate the weight functions. Especially 



the Wigner formula for rotation matrices [3i|, |3^ is crit- 
ical. In this paper we give the explicit forms of weight 
functions for j — 1/2,1 and 3/2 (two-, three- and four- 
component models, respectively) and these results im- 
ply that the weight functions Ai are generally given 
by polynomials. Through these polynomials, the initial- 
qudit dependence of the limit distribution of pseudove- 
locity is completely determined. 

This paper is organized as follows. In Sec. II, Wigner 
formula of (2j -I- l)-dimensional irreducible representation 
of the rotation group SO (3) is summarized. The fam- 
ily of quantum-walk models associated with the rotation 
matrices is defined for quantum walkers with (2j -1-1)- 
component qudit in Sec. III. In Sec. IV we also introduce 
the tensor-product models of one-dimensional quantum 
walks associated with the rotation operator R follow- 
ing the general theory of Brun et al. 18]. There we 
show the reducibility of the tensor-product models to 
our models. Sec.V is devoted to proving the general- 
ized weak limit-theorem (convergence of all moments) for 
pseudovelocities of quantum walks. There the polynomi- 
als, which give the weights of Konno's density functions 
and the point mass in the limit distributions, are listed 
for j = 1/2, 1 and 3/2, explicitly. Comparison with com- 
puter simulation with the present analytic results is given 
in Sec. VI. In this last section, we also discuss relations 
between our results and other multi-component models 
[Ulli, H m, m, Hi] and possible future problems. 



m, m' = — J, — j + 1, ■ ■ ■ , j . We can show that 

Rlil («, f3, 7) = e-'-'-r^L' (/3)e-'-^"' (6) 

with 



cos ■ 



2j+m-m' -2£ 



sm ■ 



/3 



— m 



(7) 



where 



T{j,m,m',e) 



^ Y \/ [j + m)\{i - m)\{j + m')\[j - m')\ 
^ ' {j - m' - + m - t}W.{i + m' - m)V ^ ' 

In ([7]) the summation extends over all integers of £ for 
which the arguments of the factorials are positive or null 
(0! = 1). The matrix ^ gives a (2j -I- l)-dimensional 
irreducible representation of the rotation group SO (3) 
and is called the rotation matrix. Eq.([7]) is known as the 
Wigner formula (3ll. [3^. In the present paper, when we 
write matrices and vectors whose elements are labeled 
by m,m', we will assume that the indices m and m' run 
from j to —j in steps of —1. In Appendix A, we give 

exphcit expressions of matrices = (r^^rn'^l^)) ^^"^ 

j = 1/2,1 and 3/2. 



II. WIGNER FORMULA OF 
(2j + l)-DIMENSIONAL REPRESENTATIONS OF 
ROTATION GROUP 

Any rotation in the three-dimensional real space is 
uniquely specified by three rotation angles a,/?, 7 called 
the Euler angles. In the quantum mechanics, the rotation 
with the Euler angles a,(3,j is given by an operator of 
the form (see, for instance, [3lj) 

7?(a,^,7) ^ e-'"•^3e-'^•^^e-''^^^ (3) 

where J ~ (Ji, J2, J3) is the vector operator of angular 
momentum, whose elements satisfy the su(2) Lie algebra, 

3 

[Jfc, if] i ^ £fcf„iJ„i, fc,£=l,2,3, (4) 

m— 1 

with the completely antisymmetrical tensor with three 
indices Bkim (28j . Let the ket- vector Jj, to),^ = 
0, 1/2, 1, 3/2, • • • , m — —j, — -f- 1, • • • , j, denote the nor- 
malized eigenstates of = Y^k=i ^^"^ "^3 ^^^^ that 
P\j,m) = j{j + l)\j,m), and Jslj^m) = m\j,m). (We 
set ?i = 1 in this paper.) Then, for each fixed value 
of half- integer j, a [2] + 1) x (2j + 1) unitary matrix 

R'^^\a,(3,j) = (i?w("'/3,7)) is defined with its ele- 
ments 

^w(«-/3-7) = {j,m\R{a,(3,j)\j,m'), (5) 



III. QUANTUM- WALK MODELS WITH 
(2j + l)-COMPONENT QUDITS 

Here we propose a new family of models of quantum 
walks on the one-dimensional lattice Z, in which each 
walker has a {2j + l)-component qudit ([1]). In the pre- 
vious paper : 25] , we reported the weak limit-theorem for 
the two-component model. That model is generated by 
a quantum coin represented by a matrix in SU(2) 

ue[-l,l],e,c^e[~n,7r) (9) 

and a spatial shift-operator on Z, which is represented 
by a matrix 

in the fc-space. If we compare these matrices with ^ 
with j = 1/2 and (jAip in Appendix A, we find that they 
are the special cases of i?'^/^-'(a, /3, 7); 

A = i?(i/2)(7r-6'-0,2arccos(M),-7r-6' + 0), 
S{k) = ii'(i/2)(-2fc,0,0). (11) 

From the view point of the group theory, we can give the 
following remark. In [25] we used the fact that SU(2)~ 



4 



(= the three dimensional unit sphere in R"^) and the 
quantum coin A was parameterized by three real num- 
bers, u,9,(t) (the Cayley-Klein parameters), correspond- 
ing the dimensionality 3 of the group space. On the other 
hand, we are now regarding the quantum coin ^4 as a two- 
dimensional representation of the rotation group SO (3), 
and thus the three parameters are identified with the 
Euler angles for rotations in the three-dimensional real- 
space R'^. 

This observation had led us to adopt the {2j + 
l)-dimensional representation of the rotation group, 
i?^-'^(a,/3,7), as a quantum coin to mix (2j -|- 
1) components in qudit ([T]). The spatial shift- 
matrix is given by S^J^k) = (-2fc, 0, 0) = 

We assume at the initial time t — that the walker is 
located at the origin. Then, in the fc-space, the (2j + 1)- 
component wave function of the walker at time t is given 
by 

*W(fc,t)- t = 0,l,2,..., (12) 

where 

= R^^\a-2k,f3,j). (13) 

The time evolution in the real space Z is then obtained 
by performing the Fourier transformation. 




#0')(fc,i) = ^*(^)(x,i)e-"''^, (14) 

as 

i 

¥^Xx,t+l)^ J2 i?HL'*E'(2^+2m,i),t-0,l,2,.--, 

m' — —j 

(15) 

where "im\xj t) denotes the m-th component of the (2j + 
l)-component wave function ^'(•')(x,i). 

Now the stochastic process of (2j-|-l)-component quan- 

tum walk is defined on Z as follows. Let be the 

position of the walker at time t. The probability that we 
find the walker at site x S Z at time t is given by 

Prob(Xp^ ^x) = P{x,t) = [*(^)(a;,<)]t5'(j)(a;,t), (16) 

where [\I''-')(a;, i)]^ is the hermitian conjugate of 
\l/(-')(a;,t). As shown in [lil, the r-th moment of x[^'' 
is given by 

./;j[io.,M,i.(4)'i<..,M, 

r = 0,l,2,-.-. (17) 



IV. TENSOR PRODUCT MODELS 
ASSOCIATED WITH ROTATION OPERATOR 
AND THEIR REDUCIBILITY 

A. Tensor product models 

In this section we use the notation |1) = |l/2,l/2), 
I - 1) = |l/2,-l/2) for the binary states with j = 1/2. 
Let M e {2,3,4, •••}. For an M-component variable 

m = (mi, 777.2, • ■ • , m-Af ) with 777„ G { — 1, 1}, 1 < 77 < M , 

we will write |m| = X]n=i 

Following Brun et al. [1^ we consider the 
2*^-dimensional space spanned by the bases 
{|m)}jng{_i ijAj , which are defined as tensor prod- 
ucts 

M 

|m) = n I™")- (18) 

ri=l 

Note that they are orthonormal; (m|m') = (5m, m' = 

On— 1 ^''^n'ni'n ' Let 

Rm!L'(«-/^-7) = (m|i?(a,/?,7)|m'), 

and define a 2^^ x 2*-^ matrix R[*^1 (a, /3, 7) = 

(Rm!L'("'/5'7))m,m'G{-i,i}*f- % definition of tensor 
products [Hi, 

M 

<L(«,/3,7) = WR^'X^u^i^^M) 

n=l 

M 
n=l 

where A^/'^\(3) = (r-lil/2,™72)™,™'e{-i,i} is given by 
(lAip in Appendix A. This gives a 2''^-dimensional tensor- 
product representation of the rotation group. De- 
fine the shift matrix in the fc-space by sj^^j^j/ (fc) — 
Rm!L'(-2fc,0,0) = e'^^I^IW, and set 

Y[Af](fc) ^ S[*^l(fc)RM(a,/3,7) 

= R[^^l(a-2fc,/?,7). (19) 

Following the general theory of Brun et al. [18], the wave 
function in the fc-space at time t of the one-dimensional 
quantum walk associated with the above tensor-product 
representation of SO (3) will be given by 

$M(fc,i) = (vM(fc))Vr, i = 0,l,2,..., (20) 

where the initial state is given by the 2Af-component 
qudit 

^^^^=®{q^)^ Qt.Qn^C (21) 
n— 1 ^ 
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with the normahzation Yln=ii\Qt? + \Qn?) = 1- The 
real-space wave function is then given by its Fourier 
transformation 



27r 



(22) 



Let yj*^! , t = 0, 1, 2, • • • , be the position of the walker 
at time t of this tensor-product model. The probability 
distribution function is defined by 



(23) 

The initial position of this quantum walk is the origin, 
yj*^' = 0, and the walker has the initial qudit (HI]). 



examples in the following subsection.) For (fTQ]) it implies 
that 

(V[*^l(fc))* - (K[*^1)-i 0(1/0'^.) (fc))*KM, (27) 

where V^^'^^\k) is the ^^-th copy of ([1 
Let 



and define 



(28) 



(29) 



Then it is easy to prove that the probability distribution 
function ([23]) is decomposed as 



B. Reduction of quantum- walk models 

Irreducible representations of the rotation group 

are given in the spaces spanned by \j, m), i = 

0,l/2,l,3/2,..., m = + mil. 

The two kinds of bases {|m)} and } are related 

through 



with 



j lj = l rrij 



(24) 



(25) 



where dj is the multiplicity of the (2j-|-l)-dimensional ir- 
reducible representations included in the 2^^-dimensional 
tensor-product representation and for each j an index £j 
runs from 1 to dj. Remark that J^j'^ji'^J + 1) = 2*^. 
Define the 2*f-dimensional matrix KI*-'! (K^*^' ), 

which is an orthogonal matrix; (K^*^!)^^ =*K[^1. 
Then we see 

KMR[A/l(«_2fc,A7)(KM)"' 

3 

©i?0-.d.)(c-2A:,/3,7)}, (26) 

where R^^'^^\£ = 1, 2, • • • , dj, are dj copies of the (2j + l)- 
dimensional irreducible representation of rotation group 
explained in Sec.IL That is, R[*^l(a - 2fc,/3,7) can be 
block-diagonalized into a direct sum of rotation matrices. 
Note that direct sums in ((26)) and equations below are 
taken only over i's such that k'*^' ^ 0. (See the 



pM(a;,t) = J2J2 P^''''^P^''''\^.t), (30) 

where P^^'^^^x, t) is the probability distribution function 
of our quantum-walk model introduced in Sec. Ill, whose 
initial (2j + l)-component qudit is given by 



(31) 



By this formula, the probability laws of quantum walks 
of the tensor-product models are completely determined 
by those of the models studied in the present paper. 



C. Examples 

M = 2 case We set the 2^ 



4 states 



{|(mi,TO2))}(„ii,m2)G{-i.i}2 in the order |(1,1)), 
1(1,-1)), |(-1,1)), |(-1,-1)). Then we have 
R-m' m' (", /?, 7) = e-'"(™i+"'^)~''^(™'i+™2)r|^l^, with 



-cs —cs s 



cs c 
zs 

o2 



cs —s^ c 



(32) 



r[2l = pl \ 

K m,m' J 

CS cs c 

The shift matrix is given in the fc-space by 

/ e2"= \ 

10 

1 

\ e-2ifc y 

From the above four states, {\j,m) : j = 1,0} are ob- 
tained by the highest weight construction (see, for ex- 
ample, [28| ) and we find the transformation matrix K^^^ 
as 



S[2l(fc) = R[2l(-2fc,0,0) = 



(33) 



K'^I = (^{j,m\{mi,ni2)) 



/ 1 \ 

yT72 

1 

V vX/2 J 



(34) 
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It is easy to confirm that 



K[2lr[2l(/3)(K[2l)-i 



/ 


C2 


-V2cs 


S2 








2c2-l 


-\/2cs 









%/2cs 







V 











1/ 



where r*^^^ is given by (jAip in Appendix A. This im- 
pHes K[2lvPl(fc)(K[2l)-i = y(i)(fc) © 1. This decomposi- 
tion wiU be symboHcally denoted as 



2(g)2 = 3©l. 



(36) 



= r(i)(/3)el, 



(35) 



M = 3 case We set the 2^^ = 8 states {\{mi,m2,m3))}(^jni.m2,m3)e{-i.i}^ in the order |(1,1,1)), |(1,1,— 1)), 



1(1,-1,1)), 1(1,-1,-1)), |(-1,1,1)), |(-1,1,-1)), |(-1,-1,1)), |(-1, 

g — ia(m 1+7712+7713 )— 17(771^ +7712 +'m3)^ [3] 



"'m.m' with 



-1,-1)). Then we have R^^^,ia,p.j) 



r[3] 



.[3] 





c3 


2 

— C S 


2 

— C S 


2 

CS 


2 

— CS 


2 

CS 


2 

CS 


-S3 


\ 


2 
C S 


C3 


2 

— CS 


2 

— CS 


2 

— CS 


2 

— CS 


S3 


2 

CS 




c^s 


2 


C3 


2 

— C S 


2 

— CS"^ 


S3 


2 

— C^S 


2 

CS^ 






9 


2 


2 
C S 


C3 


-s3 


2 

— CS 


2 

— CS 


2 

— CS 






9 

c s 


2 

— CS 


2 

— CS 


s3 


c3 


2 

— CS 


2 

— C S 


2 

CS 






9 

CS 


2 


-S3 


2 

— CS 


2 
C S 


C3 


2 

— CS 


2 

— C S 










C^S 


-CS2 


C^S 


-cs2 


c3 


-c\s 




v 


S3 


2 

CS 


CS^ 


C^S 


2 

CS^ 


2 
C S 


2 
C S 


C3 


/ 



(37) 



The shift matrix is given in the /c-space as 



S[3l(fc) = R[3](_2fc,0,0) = 



V 



,-ife 



-ifc 



(38) 



By the highest weight construction, we find in this case that we obtain a pair of two-dimensional subspaces in addition 
to one four-dimensional subspace in the decomposition; 



2(g)2(g)2 = 4©2®2. 



(39) 



That is, the multiplicities are 



1 and d 



1/2 



The orthogonal matrix is determined as 



K[3] 




(40) 
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We can see then 



K[3lr[3](/?)(K[3l)-i 





( c3 






-s3 














\ 






c(c2 - 2s2) 
























S(2c2 - s2) 


C(c2 - 2s2) 


















1 _ 


3 


/o 2 

Vocs 


A> 2 

V3c s 


3 






























c 


— s 
























s 


c 






























c 


— s 























s 


c 


/ 


= r 


(3/2) 




9r(i/2,2)(^)^ 















(41) 

where r(3/2)(/3) is ((X3|l and both of r(i/2'i)(/3) and r(i/2^2)(/3) are identified with ((XT|l in Appendix A. This iniphes 

K[3lv[3l(fc)(K[3l)-l = X/(3/2)(^) ^ \/(l/2,l)(^) ^ 1/(V2,2)(^)^ (42) 

I 



V. LIMIT DISTRIBUTIONS OF QUANTUM 
WALKERS 

A. Decomposition of time-evolution matrix 

A key lemma for the foUowing analysis of the quantum- 
walk models is the fact that the time-evolution matrix 
V'^^^ (fc) defined by (|13p is decomposed into the three ro- 
tation matrices 's of the form 

V^'\k) ^ i?(j')(0(fc),6'(fc),O)i?(^')(-p(fc),O,O) 

x[i?«(,/.(fc),0(fc),O)]t, (43) 

where <j}{k)^ 6{k) and p{k) are related with the Euler an- 
gles a, /3, 7 and the wave number k by 



{(a-2fc)-7} 



tan -{(a - 2k) + j} 



— tan 



p{k) 



cos 0{k), 



sm ■ 



• P^^) ■ Qll\ 

sm — — sm6'(A;). 



(44) 



We give the proof of this formula in Appendix B. 

The formula (|43p means that the time-evolution matrix 
V''^^ (k) can be diagonalized to i?*^-') {—p{k) ,0,0) by a uni- 
tary transformation given by R^^'>{<j){k),6{k),0). Indeed 
(-p(fc), 0, 0) is a diagonal matrix (-p(fc), 0, 0) = 
diagfe'^PC^), e't^-i^PC^), • • • , e-'^P^^)) and, by the unitarity 
of [i?(i)]t = (HH) is written as 

¥^^\k,t) = i?(^')(</.(fc),6i(fc),0) 



Jtjp{k) 



JtU-l)p{k) 



e'*"P«vW(fc)C^)(fc), (45) 



where Vm'' (fc) is the m-th column vector in the matrix 
i?(^)(0(fc),0(fc),O), 



vH'(fc) = 



/ i?5.£(0(fc),0(fc),o) \ 

R^'\Mk),0{k),0) 

V i?L^],„(0(fc),0(fc),o) y 



and 



C(i)(fc)^[v(i)(fc)]t0«= ^ i?(;^j„(0(fc),0(fc),O)g™,, 

m' — — J 

(46) 

where z denotes the complex conjugate of a complex 
number z. 

The expansion ([IS]) gives 



eit™pWv«(fc)c2)(fc) r 



'm=-j 



Since i?^-'^ is unitary, its column vectors make a set of 
orthonormal vectors, [vm\k)]^vl^,{k) = 6mm'- Then we 
have 



E 



|c«(fc)|V + o(r-i), 



and thus (fTT)) gives the following expression for moments 
of pseudovelocity x'f^ /t in the long-time limit (Til. [25j. 



lim 

t — >oo 



X 



U) 
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E r ^{i-mc^dHkr + \d^l{k)\'} 

m:0<m<j ^ 



(47) 



r — 1, 2, 3, • • •, where the summation is taken over m = 
1/2,3/2,- • - jj, if j is a half of odd number, and m = 
1,2,- • if j is a positive integer. Here it should be 
noted that, when j is a positive integer, m — mode 
exists, but it does not contribute to any moment of order 
r = 1, 2, 3, • • • in ([T71) . The m = mode comes from the 
fact that the walker can stay at the same position in a 
step, when (2j + 1) is odd, and its contribution to the 
limit distribution will be described by a point mass at 
the origin (see Sec.V.C). 



B. Planar orbits in parameter space and integrals 

The equations (|44|) define a one-parameter family (with 
parameter k) of transformations from the Euler angles 
(a, (3, 7) to {p, 9, (p). More explicitly, we can find the fol- 
lowing equations from (|44|) (see Appendix B), 

p(k) (3 1, 

cos ^ = cos — cos -(a + 7 — 2k), (48) 

sin^ ^ ^1 - cos2(/3/2) cos2{(a + 7 - 2fc)/2}, 

(49) 

cos(/3/2) sin{(a + 7 - 2fc)/2} 
Vl - cos2(/3/2) cos2{(a + 7 - 2fc)/2y' 

(50) 

sin(/3/2) 



cos9{k) = — 



sm0{k) 



a/I - cos2(/3/2) cos2{(a + 7 - 2fc)/2} 



0(fc) = _(a-7-2fc-^). 



(51) 
(52) 



Following the argument given in [25| , we consider a vec- 
tor p(fc) — {pi{k),p2{k),p3{k)) in the three-dimensional 
parameter space defined by 

Pi(fc) — p {k) sin 9 (k) cos (j){k), 
P2{k) = p{k) sin 9 (k) sin (t){k), 
Psik) = pik)cose{k). (53) 

Let 

ei = (— sin7, — C0S7, 0), 

62 = (sm-cos7, -sm-sm7,-cos-), 

63 — (cos — C0S7, — cos — sm7, sm — ). (54) 

Using (|48 p - (|5^ . it is easy to confirm the fact that 
p{k) _L 63 for all k G [— 7r,7r), 



which implies that p(fc) draws an orbit on a plane includ- 
ing the origin, whose normal vector is 63 in the param- 
eter space. On this orbital plane, we define the angle x 
by cosx = p(^) • Gi, where p(fc) ~ p{k)/p{k). Then we 
have the relations 



cosx 



smx = 



sin(/3/2) cos{(a + 7 - 2fc)/2} 
^1 - cos2(/3/2) cos2{(a + 7 - 2fc)/2} ' 

sin{(Q; + 7-2fc)/2} 
y/1 - cos2(/3/2) cos2{(a -(- 7 - 2fc)/2} ' 



(55) 
(56) 



Comparing (|55p with (j48|) and (|49| . the equation of the 
orbit on the plane is determined of essentially the same 
form as reported in [2^ . 



tan — — — tan ■ 



2 cosx 



(57) 



As pointed by [25|, the integral with respect to the 
wave number k in (|47p is mapped to the curvilinear in- 
tegration along the orbit with respect to the angle x 
through the relations ([55)1 and (|56p . or their inverted 
forms 

1 



cos -(a -I- 7 — 2fc) = 



cos X 



sin - [a 
2^ 



7- 2fc) 



1 



The Jacobian associated with the map k 
as 

J= — 

~ dX 

From (l49l). we have 



(58) 

= . (59) 
cos2(/3/2)sin2x 

X is obtained 



1 -cos2(/3/2)sin2x 
sin X sin(/3/2) 



sin(/3/2) 



1 -cos2(/3/2)sin^X 



(60) 



pik) 
and then 



/3 1 

2 arccos cos — cos - (a + 7 — 2fc) 



^ = -2cos-smx, 



(61) 



where the formula ((i/(ia;)arccosx — ^l/y/l — x'^ has 
been used. The long-time limit ([47|) of moments of pseu- 
dovelocity is now expressed as 



m:0<m<j 



(62) 



where 



m \ 



" rfx 



sin(/3/2) 



with (72I(x) 



_^ 27rl-cos2(/3/2)sin2x 

{|Ci;?nx)P + (-i)1(5!'I(x)P} 



2m cos — sin x 



cE(fc(x)). 



(63) 
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C. Superposition of scaled Konno's density 
functions 

In each integral Im\r), we change the variable of in- 
tegral from X to J/ by 

y = 2tocos ^ sinx. (64) 

If we assume that by this change of variable (7m "* (x) is 
replaced by Cm^ (y) , the integral is written as 

x{|c(^n2/)P + (-inc^l(y)P}- (65) 

Here fj,{x;a) is given by ([2]), which is the density 
function first introduced by Konno to describe 
the limit distributions of the two-component one- 
dimensional quantum walks in his weak limit-theorem 

i, [3. As a function of \cl^ (y)P is separated 
into an even-function part and an odd-function 

part. For positive values of m, let A4iv^\y/2'm) = 

even- function part of | cin (y) P + | c^,\j (y ) ^ and 

Ml^^^\y/2m) = odd-fmiction part of |cm''(2/)p — 

\c^-ln{y)\'^- Since ^(x;a) is an even function of x, ([SS]) 
gives 




(66) 

for n = 1, 2, 3, • • •. Then ^ implies that 




for r = 1, 2, 3, • • •, where 

M (x) = (^'^) (x) + Mi'^f (x) . (68) 



When j is a positive integer, {i.e., (2j + 1) is odd), the 
integral 




(69) 



is generally less than one, since the contribution from 
the m — mode is not included in the summation. The 
difference 

A(i) = 1 - J(J) (70) 

gives the weight of a point mass at ?/ = in the distribu- 
tion. 

The result is summarized as follows. The long- 
time limit of the pseudovelocity of (2j -I- l)-component 
quantum walk is described by the probability mea- 
sure, which consists of the summation of appropriately 
scaled Konno's density functions with weight functions 
■■™-\y / 2m) , and a point mass at the origin with 
weight A*^^^, if the number of components {2j + 1) is 
odd, that is 

^lim((-^j) = y dyyVW(y), r = 0,l,2,-.. 

(71) 

with 

^"fa)= E iM(^;c„4)A<u,.,.>(^) 

m:0<m<j ^ ' 

+l{(2,+i)isodd}A(^'^'5(y). (72) 



D. Polynomials A^'-''™'(x) representing parameter 
and initial-qudit dependence 

Using the formulas given in Appendix C and the matri- 
ces r^^^ in Appendix A, the weights M.'^^^'^^ix) and A'-') 
in the limit distribution ([7^ are explicitly determined as 
follows for i = 1/2, 1 and 3/2, < m < j. Set 

T = taii^. (73) 

For a complex number z, Re{z} denotes the real part of 
z. 

j = 1/2 case (two-component model) 



^(V2,l/2)(^)^l+_^(l/2,l/2)^^ (74) 

with 

_^(i/2a/2) ^ _{|5,/2p-k-i/2in+2rRe{gi/2g_i/2e-'^}. 

(75) 

When m[^^^''^^^^ = (resp. TW^^/^^^/^' ^ 0), the prob- 
ability density function of limit distribution v^-^/^^ (y) is 
symmetric (resp. asymmetric) 0, 



j = 1 case (three-component model) 



with 



M 



(14) 



l{\qi\^ + 2\qo\' + \q-i\'} - Re{q,q_,e-^"^], 



M'i''^ = -{|gi|2-|g_ip} + %/2TRe{(9igo + 9o«-i)e-'^}, 

M^^''^ = - 2\qo\^ + \q-if} - V2TRe{{q,qo - qoq-i)e-'^} + (1 + 2t^) Rc{q,q_,e-''-'} , 



and 



The condition that the probabiUty density function of hmit distribution v^^^y) is symmetric is given by Ai'l^ 
Generally the point mass at the origin appears with the weight A^^^ in the limit distribution. 
j = 3/2 case (four-component model) 



and 



with 



^(3/2,3/2) = A^(3/2.3/2) + ^(3/2,3/2)^ ^ M^^^^'^^'^x' + M^i^^'^/'^X^ 

^(3/2,3/2) ^ i{|53/2|2 +3|gi/2|2 + 3|9_i/2|' + |9-3/2p} " ^Re{ (93/29-1/2 + 91/29-3/2)6"''^ } , 
(3/2,3/2) ^ _^||^3/,|2 + |g^/^|2 _ |^_^^^|2 _ |^_^^^|2| _ Re{g3/29-3/2e-'"' - 91/29-1/26"'^} 



A^p/2,3/2) 



^^(3/2,3/2) 



Re| (53/291/2 + 9-i/29-3/2)e + ^Rej (93/29-1/2 - 9i/29-3/2)e 
^{|93/2|' - I91/2I' - I9-1/2I' + I9-3/2I'} - \/3TRe| (93/291/2 - 9-i/29-3/2)e"'^} 
+ ^(1 + 2r2) Re{ (93/29-1/2 + 9i/29-3/2)e"'"'}, 

-J{|93/2|' -3|9i/2|' +3I9-1/2I' - I9-3/2I'} + ^T(3 + 4r2)Re{ 93/29-3/26-''^} 

V3 



-2'^I^<5{9i/29-i/2e '''} + -^TRe|(93/29i/2 + 9-1/29-3/2 
--^(1 + 2t') 1^61(93/29-1/2 - 9i/29-3/2)e-''^}, 



)e-^} 



and with 



(3/2,1/2) _ 



|{3|93/2|' + |9i/2p + I9-1/2I' + 3I9-3/2I'} + ^Re| (93/29-1/2 + 9i/29-3/2)e ''^}, 
^(3/2,1/2) ^ ^i|3|93/2p - 5|9i/2|' + 5I9-1/2I' - 3|9-3/2p} + Re{93/29-3/2e-''^} 

^^TR-c{9i/29-i/2e-''^} + ^TRe{(93/29i/2 +9-1/29-3/2)6-''^} 
-^Re| (93/29-1/2 - 9i/29-3/2)e-''^}, 

{I93/2P - |9i/2p - I9-1/2P + I9-3/2P} + \/3t Re I (93/291/2 - 9-i/29-3/2)e-'''} 



(3/2,1/2) _ 
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M 



(3/2,1/2) 



Vs. 

2 



(1 + 2t2) Re{(q3/2g_i/2 + 91/29-3/2)6-''^}, 

^{193/21' - 3|gi/2p + 3|q_i/2p - k-3/2p} - lr{3 + ir') Re{q3/29-3/2e-''^} 
' " ^)Re|(q3/2g_i/2 -9i/29-3/2)e"'''^}- 



'3/291/2 + 'Z-1/2'7-3/2 



2t^ 



(82) 



If and only ii m[^^^'^^^^ = TW^^/^'^/'^ = and m[^^^'^^^^ = TWf /''^/'' = 0, the probability density function of limit 
distribution i/^^/^) (y) is symmetric. 



r 



These results imply that A4'^^'™^(x) are polynomials of a; 



of degree 2j and the coefhcients Ai 



that the factor i is irrelevant for limit distribution, but 
by this factor is in SU(2), see [H].) When we 



and 



but they do not on a. It should be noted that 



depend on 13 and 7 through the functions r = tan(/3/2) choose the initial qubit as Vo = (1 + i, 1 - i)/2, the hmit 

distribution of pseudovelocity is symmetric as shown by 
Figl2] (a), but when we choose *(/)o = (1 + i, 1 + i)/2, 
only changing the sign of imaginary part of the second 
component, the limit distribution becomes asymmetric as 
shown by Fig. [l](b). When (a,/3,7) = (0, -37r/2, tt), the 

former initial qubit satisfies the condition A^^^^''^^'-* = 

fl/2 1/2) 

0, but the latter does not, where M\ ' ' ' ' is given by 



their dependence on initial qudit ([T]) is complicated. In 
other words, the limit distribution of pseudovelocity of 
quantum walk is very sensitive to changes of initial qudit. 



VI. COMPARISON WITH COMPUTER 
SIMULATIONS AND CONCLUDING REMARKS 

In order to demonstrate the validity of the above re- 
sults, here we show comparison with computer simulation 
results. In the following figures, Figsl2][4l the scattering 
thin lines indicate the density of distribution oi Xt/t at 
time step t = \QQ obtained by computer simulation and 
the thick lines the probability densities of limit distribu- 
tions v^-'\y) given in the previous section. Note that if 
we integrate the density over an interval [a,b], then we 
obtain the probability that the value of Xf/t is in [a, b]. 

Two-component model 






14 




12 




10 




dens 




4 




2 


ll. 







FIG. 2: (Color online) Comparison between simulation results 
and the probability densities of limit distributions for the two- 
component model, (a) Symmetric and (b) asymmetric cases. 



The result (|74p with (1751 is completely identified with 
the previous results (ol. lid. [25|. Here we put (a,/?, 7) — 
(0, — 37r/2, tt). Then we have 

«<.« = -L(;_\), ,8„ 

which is the Hadamard matrix multiplied by i. (Remark 



(|75p . The shape of probability density function in limit 
distribution is very sensitive to changes of initial qubit. 

Three-component model 




FIG. 3: (Color online) Comparison between simulation re- 
sults and the probability densities of limit distributions for the 
three-component model, (a) Symmetric and (b) asymmetric 
cases. The perpendicular thick lines at the origin indicate 
Dirac's delta functions. 

If we set (a,/3, 7) = (0, arccos(— 1/3), tt), the three- 
component quantum coin will be 



3 



-1 -2 
-2 -1 
-2 2 



(84) 



-1 



Remark that it is similar to the Grover matrix [12| . 
but it is not the same. Figure [3] shows the compari- 
son of simulation results and limit distributions for (a) 
*<^o = (1 — 1 + i, 1 — i)/\/6, which gives symmetric 
distribution, and for (b) *-(j)Q ~ (1 — i, 1 — i, 1 — i)/\/6, 
which gives asymmetric distribution, respectively. It is 
readily checked that the case (a) satisfies the condition 
A^^^'^'' = for symmetric distribution. In the three- 



12 



component model, Dirac's delta-function-type peak at 
the origin is usually observed in simulation. 

Four-component model 





FIG. 4: (Color online) Comparison between simulation results 
and the probability densities of limit distributions for the four- 
component model, (a) Symmetric and (b) asymmetric cases. 



We set (a,/3,7) = (0, 27r/3, tt) for the four-component 
model, which corresponds to choosing the quantum coin 
as 



ij(3/2) 



( 1 3 3n/3 3^3 

3 5 \/3 -3V3 

3\/3 %/3 -5 3 

V3\/3 -3V3 3 -1 y 



(85) 



If we assume the initial qudit as 



2%/5 



/ l + 3i 



V -3-Fi 



and 




(86) 



the distributions are determined as shown in Fig. 3] (a) 
and Fig. |l](b), respectively. 

We observe oscillatory behavior of density functions of 
Xt/t in computer simulations. In general, as the time 
step t increases, the frequency of oscillation becomes 
higher. But the convergence of any moments given by 
((7T|) means that, if we smear out the oscillatory behav- 
ior, the averaged lines of density functions will be well- 
described by those of the limit distributions ([7^ . which 
is the fact demonstrated by the above figures. 

Now we discuss relations between our models and other 
multi-component models and possible future problems. 
Inui and Konno [33] and Inui et al. [34] introduced 
a one-dimensional three-component quantum walk and 
showed that a kind of localization phenomenon occurs in 
their model by calculating the long-time distribution of 
walker's position Xt- They claimed that such a localiza- 
tion phenomenon results in a point mass at the origin, 
represented by Dirac's delta function, in the limit distri- 
bution of Xt/t. Though their model associated with the 



Grover matrix does not belong to our family of models, 
the structure of probability density function in limit dis- 
tribution obtained in our three-component model [j = 1 
case) (17^ with j — m — 1 and with ([75)1 is very similar 
to the limit density function of Xt/t, given by Eq.(16) in 
[s^ ]. Then the present work suggests that such a local- 
ization phenomenon is universal for the models, in which 
there is a positive probability for walker to stay at the 
same position in each step. Further study of localization 
phenomena in quantum-walk models will be an interest- 
ing future problem. The relation between the present 
models and the tensor-product models of Brun et al. [l^ 
was already discussed in Sec. IV. Similarity between the 
density functions of walkers plotted in Fig.l in their pa- 
per and ours shown above will be explained by reducibil- 
ity of tensor-product models. As demonstrated in Sec. IV. 
C, their 2^-dimensional model will include our three- 
component model and their 2'^-dimcnsional model does 
our four-component one, and thus the density function 
shown by Fig. 1(b) (resp. Fig. 1(c)) of Brun et al. |l3] has 
the same structure with Fig[3Kb) (resp. FiglDJb)) in the 
present paper. Venegas-Andraca et al. (30| also reported 
quantum-walk models with entangle coins, where a vari- 
ety of distributions of walker's position with multi-peak 
zones are plotted in figures. They constructed quantum 
coins for multi-component qudits by also considering ten- 
sor products of the basic two-component quantum coins. 
Systematic classifications of multi-component quantum- 
walk models and patterns of their limit distributions will 
be an important future problem. In the present paper, we 
gave the explicit expressions of M^^'^"'^{x) for j — 1/2, 1 
and 3/2. Here we gave only key formulas in Appendix 
C for calculation of them, but the present result im- 
plies that Ai^^'"^''{x) are generally given by polynomi- 
als. We hope that further study of is a promising 
subject in the field of quantum walks. As reported in 
[ill . [itI . [19I . [2ll . [sol . [35} . multi-component models have 
been studied to simulate quantum walks on the plane, in 
the higher-dimensional lattices Z'^, or on general graphs. 
The present work suggests that the group-theoretical in- 
vestigation will be useful to perform systematic study of 
such extended models of quantum walks and quantum 
processes. 
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APPENDIX A: THE MATRICES A^\l3) FOR j = 1/2, 1, 3/2 

The explicit expressions of A^\j3) are readily derived from the Wigner formula ([7]) as follows for j — 1/2,1,3/2. 
Let c = cos(/3/2), s = sin(/3/2). Then 



r(i/2)(/3) 

^(3/2) 



c — s 
s c 



V2cs 2c2 - 1 -V2c 



V2i 



cs 



( <? 



(Al) 
(A2) 



(A3) 



APPENDIX B: PROOF OF EQUATIONS (|43]) WITH (|44l) 

By p^. it is enough to prove the equality 

(a, /3, 7) = (0, 61, 0)i?(^') 0, 0) [R^^'^ (0, 0, 0)]^ 

with relations 

1, .,7r 1 

— (a — 7) = (pH — , tan — (a + 7) = — tan — cos d, sm — — sm — sm fy. 
The equality (jBl[) is a matrix representation of the equality which will hold among the rotations 

/3, 7) = e, 0)R{-p, 0, 0) [i?(</), 0, 0)]t . 



(Bl) 



(B2) 



(B3) 



Since rotation matrices defined by ([5]) with j = 1/2, 1, 3/2, • • • are faithful representations of rotations, it may be 
enough to prove the equality (IBip for the simplest case, j — 1/2, since it implies (|B3p . By direct calculation, we have 



0^ 0)i?(i/2) 0, 0)[i?(i/2)(0, 0, 0)]t 
e-''^/2cos(6'/2) -e-''^/2gij^(5,/2) \ / 
ei0/2gin(6i/2) e'-^/^ cos(6i/2) j V ^-'p/^ 

cos(p/2) + i sin(p/2) cos i sin(p/2) sin Se"'"^ 
i sin(p/2) sin 0e"^ cos(p/2) — i sin(p/2) cos 



e''^/2cos(6l/2) e-'"f'/^sm{0/2)\ 
-e''^/2sin(6l/2) e-'-^/^ cos(6l/2) j 



(B4) 



Comparing it with 



we have the equations 



fi>(V2)r^ « f e-'("+^)/2cos(/?/2) sin(/3/2) 

^ ia,/^,7j-^ gi(a-7)/2 sin(/3/2) e'("+^)/2 cos(/3/2) 



cos - = cos — COS - (a + 7), sm - cos = ~ cos — sm - (a + 7), 
2 2 2 2 2 2 



sin — sin sin ( 
2 



sm — cos — [a — 7], sm — smt' cos (p — sm — sm — — 7J, 

Zj Zi Zi z 



(B5) 



from which (IB2D is derived 
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APPENDIX C: FORMULAS 



Inserting ([6]) with ^ into (|46|) and taking square of the complex variable, we have the expression 

mi = ~jm2 = -j ii (.2 

/ /i/i \ \ 4j+mi+m2-2m-2^i— 2^2 / L\ \ 2f i+2£2+2m-mi -m2 

X cos-V^ ™4e ■ (CI) 



From ([50)) - (|52p . the following relations are derived 



cos 9{k) = — cos ^ sin x, 

sin6'(fc)e''^('') = ^sin^ sinx - icosx ) e"'''. (C2) 



Then, through the change of variable (|64p . we have 

2 V 2mJ ' 



„2 ^(fc) _ y 



sm = - 1 H , 

2 2 V 2m/ ' 



sin^cos^e'^t^) = J tan ^ - i./l - m%ec^ | U"- (C3) 

2 2 2 1 2m 2 V V2m/ 2 J ^ ' 

By (|C3p we can perform the transformations \dii\k)\^ ^ \Cn? {x)\^ ^ |cm and 7W(^''"')(y/2m)'s are obtained. 

In order to calculate (|78p . we have used the following integral formulas, 

2cos(,3/2) ^ 
.2cos(/3/2) ^ ^ 

Jo {1 - (2//2)2}^os2(/3/2) - (y/2F'^^ " sin(/3/2) ' ^^^^ 
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